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ABSTRACT 

Gravitational waves (GWs) from the inspiral of a neutron star (NS) or stellar- 
mass black hole (BH) into an intermediate-mass black hole (IMBH) with mass 
M ~ SOMq to 35OM0 may be detectable by the planned advanced generation 
of ground-based GW interferometers. Such intermediate mass ratio inspirals 
(IMRIs) are most likely to be found in globular clusters. We analyze four pos- 
sible IMRI formation mechanisms: (1) hardening of an NS-IMBH or BH-IMBH 
binary via three-body interactions, (2) hardening via Kozai resonance in a hier- 
archical triple system, (3) direct capture, and (4) inspiral of a CO from a tidally 
captured main-sequence star; we also discuss tidal effects when the inspiraling 
object is an NS. For each mechanism we predict the typical eccentricities of the 
resulting IMRIs. We find that IMRIs will have largely circularized by the time 
they enter the sensitivity band of ground-based detectors. Hardening of a binary 
via three-body interactions, which is likely to be the dominant mechanism for 
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IMRI formation, yields eccentricities under 10~^ when the GW frequency reaches 
10 Hz. Even among IMRIs formed via direct captures, which can have the high- 
est eccentricities, around 90% will circularize to eccentricities under 0.1 before 
the GW frequency reaches 10 Hz. We estimate the rate of IMRI coalescences 
in globular clusters and the sensitivity of a network of three Advanced LIGO 
detectors to the resulting GWs. We show that this detector network may see 
up to tens of IMRIs per year, although rates of one to a few per year may be 
more plausible. We also estimate the loss in signal-to-noise ratio that will re- 
sult from using circular IMRI templates for data analysis and find that, for the 
eccentricities we expect, this loss is negligible. 

Subject headings: black hole physics — globular clusters: general — gravitational 
waves 



Introduction 



Observational evidence from cluster dynamics and from ultra-luminous X-ray sources 
suggests that there may exist a popula tion of intermediate-raass black hol es (IMBHs) with 
masses in the M ~ 10^ - 10"^ Mq range jMiller fc Colbertlbooi iTrentilboOGf ). Numerical sim- 
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), through a combination of emission of grav- 



itational radiation, binary exchange processes, and secular evolution of hierarchical triple 
systems. Gravitational waves (GWs) will be generated during the intermediate-mass-ratio 
inspiral (IMRI) of a stellar- mass object (black hole [BH] or neutron star [NS], since a 
white dwarf or a main-sequence star would be tidally disrupted) into an IMBH. For IMBH 
mass < 350 Mq, these waves are potentially detectable with the planned advanced gen- 
erati on of ground-based GW interferomet ers: Advanced LIGO and its international part- 
ners tearish fc Weis]ll999l : IPritschellboOsh . 



IMRIs will be important as probes of strong gravity and cluster dynamics due to their 
mass range and dynamical histories. For example, from Advanced LIGO IMRI data it may be 
possible to measure the quadrupole moment, Q, of an IMBH to an accura cy of AQ ~ Qkbtt, 
where Qxerr is the quadrupole moment of a Kerr BH (IBrown et al.l 120071 ). This is sufficient 
to distinguish a BH from a boson star, for which the quadrupole moment can be many times 
the Kerr value. In addition, since the formation of IMBHs in clusters seems to require short 
mass segregation timescales (see §[2]), detection of IMBH mergers and their associated masses 
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will yield information about young dense clusters and their evolution. 

In this paper we discuss the astrophysical and data analysis aspects of IMRIs. In § [2] 
we provide the astrophysical setting for IMRIs and describe the formation mechanisms. We 
estimate the typical eccentricities resulting from different capture mechanisms and find that 
inspirals will largely circularize by the time the GW frequency reaches the Advanced LIGO 
band (/gw ^ 10 Hz). We show, in particular, that three-body hardening, which is likely to 
be the dominant IMRI formation mechanism, will result in binary eccentricities e < 10~^ in 
the Advanced LIGO band. Even direct capture, which is the most likely mechanism to yield 
high eccentricities, leads to ~ 90% of IMRIs with e < 0.1 in the Advanced LIGO band. In §[3] 
we estimate an upper limit on the rate of IMRIs detectable by Advanced LIGO of up to 10 
events per year. A more sophisticated, but model-dependent, rate estimate ranges from one 
event per 3 years for NS IMRIs to 10 events per year for 10 Mq BH IMRIs. The event rate 
can be enhanced by a factor of ~ 3.5 by optimizing Advanced LIGO for detections at low 
frequenciesll] Searches for IMRIs in Advanced LIGO data will likely use matched filtering 
techniques, for which accurate waveform templates are required. In § HI we estimate that 
there will be a negligible loss in signal-to-noise ratio (S/N) if circular templates are used to 
search for IMRIs with the expected eccentricities in Advanced LIGO data. 



2. Astrophysical Setting, Capture Mechanisms, and Typical Eccentricities 

An IMBH cannot result from the evolution of a solitary star in the current universe be- 
cause even a star of initial mass ~ 10^ Mq will be reduced well below this mas s by wi nds and 



pulsational instabilities driven by metal-line opacities (cf. iFryer fc Kalogeral (120011 ). Fig. 7 



and associated discuss ion). Some IMBHs might be formed from the first, metal-free, stars 



(IMadau fc Reed 1200 ll ). but these IMBHs are unlikely to participate in multiple mergers 



with COs. Instead, we focus on the proposal that IMBHs can be produced in the current 



stars segregate to the center in less than their ~ 2 x 10^ yr lifetimes (Ebisuzaki et al. 2001; 
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, then stellar mergers can overcome mass losses and the collision 



product can reach hundreds to thousands of solar masses, presumably evolving into an IMBH. 
When supernovae start to occur, the resulting mass loss leads to an expansion of the 



'^We used the Advanced LIGO Bench code to perform this optimization: 
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cluster, which thus transitions into a more colhsionless stage of existence. From this point 
on, COs can be captured by the IMBH and generate observable GWs as they inspiral under 
radiation reaction and eventually merge with the IMBH. 

Early in the history of the globular cluster the inspiraling objects in IMRIs are likely to 
be m ~ 10 M(7) BHs, which may form a dense subcluster composed pu rely of BHs around the 
IMBH Jo'Learv et al.lbood : lO'Learv. O'Shaughnessv. fc Rasiol 120071 ) . Late in the cluster's 
history, once the BH central subcluster has largely evaporated, NSs will likely replace the 
larger BHs as the inspiraling objects. 

There are several ways in which stellar-mass COs can be captured by an IMBH. Most 
mechanisms of capture involve binaries, because the cross section of a binary is orders of 
magnitude larger than that of a single CO. 

Extensive numerical studies of binary-single interactions demonstrate that hard binaries 
(defined, e.g., so that the total energ y of the binar y-single system is negative) tend to be 
tightened by three-body interactions (jHeggie Ill975l ). These studies also show that massive 
objects such as stellar-mass BHs and IMBHs tend to swap into binaries. The most likely 
capture mechanism involves the formation of a CO-IMBH binary, which is subsequently 
hardened by repeated three-body interactions until radiation reaction becomes significant 
and the binary coalesces. 

Hardening can also occur via binary-binary interactions; unlike binary-single interac- 
tions, these can result in a stable hierarchical triple. Some frac tion of these end up in 
orientations favorable for the secular Kozai resonance (jKozailll962l ). in which the inner bi- 
nary (which contains the IMBH) periodically increases and decreases its eccentricity while 
keeping its semimajor axis constant. The periapsis distance can theref ore be quite low in 
parts of the cycle and can lead to coalescence without Newtonian recoil (IMiller fc Hamilton 
2002bl : IWenll2002l ) . although recoil from gravitational radiation will still occur (see § l3.2p . The 
importance of the Kozai resonance depends on the frequency of binary-binary interactions, 
which is unknown at present. 

In addition to these mechanisms, which usually require multiple interactions to lead to 
merger, a hyperbolic encounter at a close enough periapsis can produce direct capture via 
emission of gravitational radiation. Assuming that the IMBH does not have a significant 
radius of influence, the effective cross section for radiative captu re of the CO by an IMBH 
is proportional to M^^/^, where M is the mass of the IMBH ( Ouinlan fc Shapiro! 119871 ). 
For two-body encounters this process is likely to be important only for masses high enough 
(> 10^ Mq) that the frequency of the GWs throughout the subsequent inspiral will be below 
the sensitivity range of ground-based detectors. However, direct capture during a three-body 
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interaction could be significant (iGiiltekin. Miller. &: Hamiltonll2006l ). 



Finally, an IMBH could tidally capture a main-sequence star. If the captured star evolves 
to a CO while in orbit around the IMBH, the remnant could remain bound to the IMBH and 
ultimately spiral in via GW emission. This scenario has been su ggested as a possible explana- 



tion for the observed population of ultra luminous X-ray sources ( jHopman. Portegies Zwart. fc Alexander 



2004 : iHopman fc Porteries Zwartll2005[ ). 



Additionally, orbital energy may couple to vibrational normal modes of the inspiraling 
object in the case when the inspiraling object is an NS. In principle, the energy loss to tidal 
heating of an NS could change the inspiral trajectory, or even disrupt the NS. 

The IMRI enters the Advanced LIGO band when 

= = i >10Hz. (1) 

i.e., when the periapsis is Vp ^ IQGM/c^ = IQOOGMq/c^ for M = IOOMq. The frequency 
of the dominant quadrupolar {k = 2) harmonic in the GWs emitted at the innermost stable 
circular orbit is 

for inspirals into non-spinning BHs. 

Below, we discuss the eccentricity of an IMRI at the time its GW frequency enters the 
Advanced LIGO band for each of the mechanisms mentioned above: (1) formation of a CO- 
IMBH binary and subsequent hardening via three-body interactions, (2) Kozai resonance of 
a hierarchical triple system, (3) direct capture when a solitary CO passes close to the IMBH, 
and (4) tidal capture of a main-sequence star that subsequently evolves to leave a CO. We 
also consider the impact of (5) tidal interactions with an inspiraling NS. 



2.1. Hardening of a CO-IMBH Binary via Three-Body Interactions 

This mechanism proceeds as follows. The IMBH rapidly swaps into a binary because 
it is far more massive than any other object in the globular cluster. Advanced LIGO can 
detect IMRIs at redshifts z < 0.2 (§ 13.11) . so it will predominantly observe globular clusters 
late in their history. On a timescale that is short relative to the merger time, an NS or a 
BH will encounter the binary containing the IMBH and will exchange for the companion 
in this binary, since stellar remnants are the most massive objects in the late cluster other 
than the IMBH itself. From that point on, stars of all types (although biased towards the 
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heavy ones that segregate towards the center) engage in three-body interactions. Numerical 
simulations show that interactions tend to tighten a binary if it is hard. This can be under- 
stood heuristically for three equal-mass objects by noting that ejection will tend to occur at 
roughly the binary orbital speed; hence, if this is greater than the initial encounter speed 
at infinity the binary loses energy. Binary tightening proceeds until the binary can merge 
through radiation reaction from the emission of GWs. 



For the unequal mass binaries we consider here, simulations by lQuinlaru (119961 ) show that 
a single interaction of a field star of mass m^, with a binary of total mass M will on average 
change the binary energy by a fractional amount AE/E = 0{m^/M), roughly independent 
of the component masses of the binary. More precisely, approximately (2tt / 22) M/ m^ inter- 



actio ns are required to reduce the semimajor axis of a hard binary by one e-folding (iQuinlan 



1996|). 

The rate at which objects interact with the IMBH binary is 

N = nqv, (3) 

where n is the number density, v is the relative speed, and <^ is the gravitationally focused 
cross section <^ = 7ra(2GM/w^) for an interloper to approach within the binary's semimajor 
axis a of the binary. Since this rate is proportional to a, the total time for the binary to 
harden until the semimajor axis equals a is dominated by the last e-folding time: 

27r M 1 „ /10^3 cm\ 

Tharden ~ 7^ — - ~ 2 X loM yr, 4 

22 m^N \ a J 

where we set m* = 0.5 M©, v = 10^ cm s~^, and n = 10^'^ pc~ ^ (the number density of all 



stars in a core-collapsed globular cluster; iPryor fc MeylanI 119931 ). 



For a sufficiently massive BH, a cusp is formed and the interactions are no longer 
described by individual binary-single encounters. We can estimate roughly the mass above 
which this occurs. Consider a core of number density ricore and velocity dispersion a. For an 
IMBH of mass M, the radius of infl uence (inside of which t he IMBH dominates the potential) 
is Tinfl = GM/a"^. For a true cusp, Bahcall fc Woli (jl976 ) showed that the number density 



would go as n(r) = ncore('"/''"infi) The total number of objects inside ri^n is then 

N{r < rir^fi) = 47rr^ncore (r/rinfl) "^/"^cir = (167r/5)rfnflracore • (5) 

^0 

Scaling to canonical values, this gives 

N{r < Tinfl) ^ 0.3(M/100 MQ)^{a/10 km s-^)~^(ncore/10^-^ pc^^) (6) 
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Therefore, in the mass range most relevant to Advanced LIGO, it is unlikely that there will be 
a significant cusp, hence our treatment of isolated binary-single interactions is reasonable. 
For more massive BHs a cusp might form, although we note that for M < 1000 Mq the 
typical distance wandered by the IMBH in the core is larger than the radius of infiuence; 
hence, cusp formation could be made more difficult. This is, however, worth further study. 

The gravitational radiation merger timescale for an intermediate-mass-ratio binary of 
semimajor axis a, eccentri city e, reduc ed mass fi ^ m, and total mass approximately equal 



to the IMBH mass M, is flPetersI 119641 1 



M^m V 10^'^ cm ' 



Simulations and phase-space arguments show that three-body inter actions cause the ec- 
centricity of the binary to sample a thermal distribution P{e)de = 2ede (iHeggie 1 ll975h in the 



Newtonian realm where gravitational radiation is not significant. If an interaction leaves the 



binary with a high eccentricity, however, it is more likely to merge. iGiiltekin. Miller, fc Hamilton 



(120061 ) examined the eccentricity of the binary after its final three-body encounter and found 
a typical value of e ~ 0.98 due to this bias. Taking this as the typical value for eccentricity, 
we find 

T ^ 10« (^] /100M,y f^^V (8) 

In fact, there is a distribution of eccentricities after the last encounter, rather than a single 
eccentricity value of 0.98. However, Monte Carlo simulations, which are described below, 
confirm the typical merger times and final eccentricities computed here analytically by as- 
suming the final-encounter eccentricity of 0.98. Moreover, simulations indicate that the 
fraction of direct plunges from highly radial orbits must be extremely small, because they 
would require improbably small periapsis separations. 

The IMRI rate will be maximized when the total merger time, T = Tharden + Emerge, 
is minimized. Minimizing T with respect to a, we find that the total merger time is T ^ 
3 X 10^ yr for the inspiral of an m = 1.4 Mq NS into an M = 100 Mq IMBH, yielding an 
IMRI rate per globular cluster of ~ 3 x 10"^ yr~^. If the inspiraling object is an m = 10 Mq 
BH, and the IMBH mass is again M = 100 Mq, then the total merger time is T 2 x 10^ yr, 
and the IMRI rate per globular cluster is ~ 5 x 10~^ yr"^- 

These numbers are close to the answers obtained with Monte Carlo simulations using the 
same procedure as in Giiltekin, Miller, & Hamilton (2006). We find from these simulations 
that the total time to merger averages 5 x 10^ yr for 1.4 Mq NSs and 3 x 10^ yr for 10 Mq 
BHs, interacting with field stars of mass 0.5 Mq and an IMBH of mass 100 Mq. We also find 
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that, as we assumed, once a CO is in a binary with an IMBH it stays there; only a fraction 
^ 6 X 10"'* of encounters swapped out an NS, and only 1 of the 5 x 10^ encounters swapped 
out a BH. Therefore, as we indicated, the object that eventually merges with the IMBH is 
highly likely to be a CO. 

This mechanism requires the cluster to be in a core-collapsed state and for this state 
to persist for a time S> 2 x lO^yr. Core collapse is expected to persist in the absence of 
significant heating, as will be the case for clusters with IMBHs in the mass range of interest, 
so 2 X lO^yr should be easily achievable. About 20% of clusters currently are in a state of 
core collapse, so this state can indeed be sustained for times of order a Hubble time, or much 
longer than 2 x lO^yr. We consider only these core-collapsed systems as likely hosts of IMRIs 
when computing event rates below. 

Radiation reaction from GW emission dominates the evolution once the GW merger 
timescale Tmerge [Eq- dH])] is shorter than the average time between three-body encounters, 
1/7V, defined by Eq. [31 For the NS-IMBH system (m = 1.4 Mq, M = 100 M©), this 
occurs when the semimajor axis takes the value a ^ 5 x 10^^ cm. As discussed earlier, the 
eccentricity at this time is e ~ 0.98, and hence the periapsis is ^ 10^^ cm ^ 7000 GM/c^. 
For the BH-IMBH system (m = 10 Mq, M = 100 Mq), radiation reaction dominates for 
a < 8 X 10^^ cm, corresponding to a periapsis of ^ 1.6 x 10^^ cm ^ 10000 GM/c^. 



Keplerian orbits evolving under radiation reaction satisfy [see Eq. (5.11) of lPeterd (Il964j )] 

rpe-i2/i9(l + e) [1 + (121/304)e2] = constant, (9) 

from which we can obtain the eccentricity at a particular frequency, given the initial values 
of periapsis and eccentricity. We find that for this capture mechanism, the eccentricity when 
the source enters the Advanced LIGO band (/gw = 10 Hz) is very small: e < 3 x 10^^ for 
the NS-IMBH system and e < 2 x 10~^ for the BH-IMBH system. The orbit will thus have 
circulariz ed by the time the IMRI is in the Ady anced LIGO band. This is consistent with the 



results of iGiiltekin. Miller, fc Hamilton! (120061 ). who also found that IMRI binaries formed 



through this channel would circularize before they entered the Advanced LIGO band. 



2.2. Kozai Resonance 



A stable hierarchical triple system could experience Kozai resonance th at would drive 
the eccentricity of the inner binary to a value close to unity (IKozail Il962l ) , leading to a 
small periapsis sepa ration and binary tighteni n g and even tual merger through gravitational 
radiation reaction (IMiller &: Hamilton! l2002bl : IWenI |2002| ). Some simulations (e.g., those 
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of lO'Leary et al.l (120061 )) suggest that the four-body (binary-binary) interactions that are 
required to place the binary on the Kozai merger track constitute only a small fraction of 
the total number of merger events in the cluster. If so, four-body interactions play a minor 
role in IMRI for mation. These simula tions may not consider all possibilities, however. In 
particular, in the lO'Leary et al.l (120061 ) model, binaries are only destroyed (through mergers, 
or by being kicked out of the subcluster). Therefore, the binary fraction decreases with time, 
meaning that binary-binary interactions are uncommon late in the cluster's history. There 
may be a way to replenish the number of BHs in binaries, however. Approximately 5% - 



20% of normal stars in globulars are in binaries (iRubenstein fc Bailynlll997l : iBellazzini et al. 



20021 ) (this fraction is closer to 50%- 70% in the field, but in globulars the wide binaries are 
disrupted). If such a binary goes through the BH subcluster, a BH could swap in, so that 
even if no BHs were originally in binaries, the binary fraction could increase. 

Although computing the relative contribution of Kozai resonance mergers to the total 
number of IMRIs requires more detailed modeling of the cluster dynamics , it is possi ble to 
estimate the largest eccentricity that could result from this mechanism (see IWeru ( 12002| ) for a 
more detailed discussion in the context of stellar-mass BHs). For this calculation, we assume 
that the Kozai resonance drives the binary to a sufficiently high eccentricity to allow merger 
via radiation reaction within one Kozai cycle. In reality, the semimajor axis and eccentricity 
would evolve gradually over multiple Kozai cycles, leading to larger typical periapses and 
smaller eccentricities, so our assumption will overestimate the typical eccentricities of IMRIs 
in the Advanced LIGO band. 

We assume that the eccent ricity is near i t s max imum for a fraction 0.01 of the total 
Kozai cycle (based on Fig. 1 of llnnanen et al.l (119971 )) and compare this time with the ra- 
diation reaction timescale. If the radiation reaction merger time is much longer than the 
time near maximum eccentricity, we assume that gravitational radiation is insignificant. If 
instead the timescale for Kozai resonance to drive the eccentricity to some value e ~ 1 is 
much larger than the timescale for radiation reaction to circularize the orbit down from e, 
then the eccentricity will never reach e in practice, even though e may be less than the 
maximum possible eccentricity for the given configuration (see below). Therefore, the max- 
imum eccentricity reachable when including gravitational radiation is given approximately 
by the condition that the radiation reaction timescale is equal to the time near that high 
eccentricity. 



The timescale for the Kozai cycle is given by, e.g., Eq. (4) of iMiller fc HamiltonI (j2002bl ): 



TKc 



few X 



m2a\ 



1/2 



Cm-; 



1/2 



3x 



Ml 



100 Mq 



1/2 



Mr. 







m2 



ai 



10^3 cm 



3/2 



(10) 
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where, in the notation of iMiller &: HamiltonI (j2002bl ). Mi is the total mass (approximately 
equal to the mass of the IMBH), mi is the mass of the inner companion, 1712 is the mass of 
the outer companion, ai is the semimajor axis of the inner binary, and 62 is the semiminor 
axis of the outer binary. Setting the timescale for merger by gravitational radiation, given 
in Eq. ([7]), equal to tgr = 0.01 XKozai yields 



ai \5/2 



1013 cm 



) 



5^/2 ~ 3 X 10-15 



Ml 



100 Mq 



5/2 



nil 
m2 



(11) 



where e = 1 



Relativistic precession constrains the maximal eccentricity, or mini mal e, that can be 
achieve d in a Kozai cycle. That minimal e is given by Eqs. (6) and (8) of iMiller fc Hamilton 
f |2002bh as: 



9 V m2afc^ 



1.6 X 10" 



1712 
Mq 



Ml 



100 M 



o 



ai 



1013 



cm 



ai 



(12) 



In order to compute the maximal plausible eccentricity at /gw = 10 Hz, we need to 
estimate the minimal periapsis radius at the peak of the Kozai cycle, when radiation reaction 
becomes dominant, since eccentricity will be close to unity there [cf. Eq. ([9])]. That is, we 
must minimize Vp = ai(l — e) ~ aie/2. This minimum value is found by solving Eqs. ( fTTI) 
and 1^. We find 



aie 



1013 



> 1.8 X 10- 



cm 



1712 

Mo 



-4/9 



Ml 



100 Mq 



13/9 



mi\ 

7712/ 



2/9 



(13) 



Stability requires that the semiminor axis of the outer binary is at least a few times 
greater than the semimajor axis of the inner binary, so we set 62/^1 = 5. We again assume 
Ml = 100 Mq, mi = 1.4 Mq, and m2 = Mq (al though this choice violate s the restricted 
three-body assumption under which Eq. (8) of flMiller fc HamiltonI l2002bl ) was derived). 
These parameter values predict a minimal Vp > 170GM/c^ at the time when radiation 
reaction takes over; hence, according to Eq. the maximal eccentricity of IMRIs formed 
via the Kozai resonance mechanism in the Advanced LIGO band is e ~ 0.01. 



2.3. Direct Captures 



If we assume that the IMBH is wandering in the stellar cluster, the effective cross 
section for direct captures via two-body relaxation (GW emission) is proportional to the 
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(12/7) power of the total mass (iQuinlan &: Shapiro! 119871 ). so an IMBH has a relatively 



small capture cross section, making this capture mechanism relatively unlikely. If we instead 
assume that the M — a relation holds for globular clusters, which is equivalent to saying 
that the IMBH dominates the dynamics in the cente r of the cluster, the capture rate would 



ap 

increase towards smaller IMBH masses, like M~^^^ ( Hopman fc Alexander] l2005l ). and this 
channel would contribute significantly to the total rate. However, as discussed in § 12. H 
the IMBHs of interest for Advanced LIGO, with M ~ lOOM©, have a very small radius of 
influence and so they will not have a significant influence on the dynamics in the cluster 
center. The direct capture mechanism, in any case, can yield higher eccentricities than 
scenarios involving binaries. 

The critical periapsis separation for the direct capture of a CO of mass m, moving at 



i nfinit y with velocity f , by an IMBH of mass M ^ m is [e.g. Eq. (11) of IQuinlan fc Shapiro 



fll989h ]: 



C^"c2 /m\2/7/ V \-4/7 

jp » Q.K,n ( — \ { \ 



GM -^^HmJ llO^^j ■ 

If M = 100 Mq, m = 1.4 Mq, and v = 10® cm s~^, direct capture is possible at a maximum 
periapsis of r^^^c^ / [GM) ^ 280; if m = 10 Mq and M and v are the same as above, 
the maximum periapsis is r^^^^c^ / {GM) ^ 500. For such small periapses, gravitational 
focusing implies oc 6^, where b is the impact parameter. Hence, the probability distribution 
P{b) oc 6 in impact parameter corresponds to a uniform distribution in periapsis at capture, 
P{rp) = constant. 

In Figure [T], we plot the eccentricity of an IMRI at the frequency at which it enters 
the Advanced LIGO band as a function of the initial periapsis at capture, following Eq. ([9]). 
The initial eccentricity at capture can be computed from the energy lost during the first 
pass; however, the exact value does not significantly affect the eccentricity at /gw = 10 
Hz, so we set the eccentricity at capture to be e = 1. The initial periapsis is uniformly 
distributed between r™^"^ = AGM/ (orbits with periapsis under 4GM/ will plunge rather 
than inspiral) and r^'^^. Therefore, to determine the total fraction of directly captured IMRIs 
that circularize to a given level e < Ccutofr by the time they are in the detector band, it is 
sufficient to find the fraction of the interval [r™™, r™*^] for which the line in Fig. [1] stays 
below Ccutoff- 

Thus, for the chosen IMBH mass of M = 100 Mq, if the CO is an m = 1.4 Mq NS, 
86% of all directly captured IMRIs will be circularized to e < 0.1 by the time they are in the 
Advanced LIGO band. If the CO is an m = 10 Mq BH, 92% of all directly captured IMRIs 
will be circularized to e < 0.1 and 67% will be circularized to e < 0.01 by the time they are 
in the detector band. 
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2.4. Tidal Capture of a Main-Sequence Star 



It has been suggested that ultraluminous X-ray (ULX) sources are systems in which 
a main-sequence sta r that has been tidally captured is transferrin g mass to an IMBH via 
Roche lobe overflow (IHopman. Portegies Zwart. fc Alexanderll2004l ). In such a system, after 
the star reaches the end of its main-seq uence lifetime and undergoes a su pernova, it may 
leave a CO on an orbit about the IMBH (IHopman fc Portegies Zwartll2005l ) and this object 
may then spiral into the IMBH via GW emission. Although work on this problem has 
focused on sources that might be detected by LISA, results have also been presented for the 
~ IOOMq IMBHs that we consider here. For M ~ IOOMq, only l%-2% of systems leave 
a CO that inspirals into the IMBH wit hin a Hubbl e time, and these remnants a re alw ays 
NSs ( iHopman fc Portegies ZwartI l2005l ). Following iHopman fc Portegies ZwartI (120051 ) we 
can estimate the rate of these events by assuming that there is ~ 1 ULX source in each 
galaxy. The ULX phase lasts approximately the main-sequence lifetime of the captured star, 
which is ~ 10^ yr, so we estimate that the capture rate is 10^'' yr^^. Multiplying by the 
fraction of events that successfully inspiral, we estimate a rate of 1-2 x 10~^ IMRIs per galaxy 
per year. There are typically ~ 100 globular clusters per galaxy, so the rate per globular 
cluster is ~ 10^^^ yr~^, which is considerably smaller than the binary hardening rate. Thus, 
while this channel could lead to some IMRIs detectable by Advanced LIGO, the rate is likely 
to be signiflcantly lower than the binary hardening channel. 

An NS cap tured via this mechanism would b egin to inspiral into the IMBH with eccen- 
tricity e < 0.9 (IHopman &: Portegies ZwartI l2005l ) and periapsis approximately equal to the 
tidal radius, {M/M^Y^^R^, where M^. > lOM© and i?* are, respectively, the mass and radius 
of the main-sequence star. Assuming, conservatively, i?* > lO^km, this capture periapsis is 
typically > 500{GM/c^). For an M = lOOM© IMBH, equation © predicts e ^ 0.002 when 
the source enters the Advanced LIGO band. In practice, the eccentricity is likely to be even 
smaller. It is thus quite clear that this capture mechanism also produces sources that are 
essentially circular when they enter the Advanced LIGO band. 



2.5. Tidal Effects 

If the inspiraling object is an NS, tides may be signiflcantly excited as it passes the 
central IMBH. If suflicient energy goes into tidal heating, the NS could be disrupted. Prior 
to disruption the orbital inspiral will be modified as orbital energy and angular momentum 
are lost into tidal heating. Tidal interactions are not important for the IMRI events we are 
considering, however, as we demonstrate below. 
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2.5.1. Tidal Disruption 

A star will be tidally disrupted by a BH when the gravitational tidal force acting over 
the star due to the BH exceeds the self-gravity of the star. Assuming a Newtonian potential, 
this leads to the usual tidal disruption radius 

f^V. 41.5k. (15) 



in which i?td is the radius at which tidal disruption occurs, m and are the mass and radius 
of the star, respectively, and M is the mass of the BH. The gravitational field outside a Kerr 
BH is not Newtonian, but (fTS!) still provides a reasonable estimate of the tidal disruption 
radius. Comparing this to the Schwarzschild radius of a lOOM© BH, = 2GM/ (? = 300km, 
suggests that, even when relativistic effects and BH spin are included, tidal disruption could 
only occur very close to the central BH. Earlier in this section we showed that the orbits 
of IMRI objects are effectively circular by the time the CO gets close to the IMBH. The 
tidal effects for stars on circular orbits are most extreme for prograde equatorial orbits, since 
these come closest to the central body. Thus, we use results for prograde, equatorial circular 
orbits for a more accurate calculation of tidal disruption. 



Vallisneril (120001 ) analyzed NS disruption using the correct tidal field for objects in 
prograde, circular, equatorial orbits around a Kerr BH and found that the GW frequency at 
which tidal disruption occurred, ftd, satisfied the relationship 



3.25km (m/1.4M0)5 (M/SOMq)^ (GM/td/c^^)""'^^ GMft^/c^ < 0.045 
1.55km (m/1.4M0)^ {M/50Mq)^ (GM/td/c^)""'^^ GMftd/c^ > 0.045 



An inspiraling object plunges into the BH when it reaches the innermos t stable prograde 



circular orbit (ISCO). This has radius (IBardeen. Press, fc Teukolskylll972l ) 

c^R- 



CM 



3 + y^Sx' + -J{3-Z){3 + Z + 2v/3x2 + Z2), 



where Z = 1 + 



(1 + X)U(1-X)^ (17) 



where x = Si/M^ is the dimensionless spin parameter of the BH. 

The condition that the star is not disrupted before plunge sets a maximum radius for the 
NS. If we require the tidal disruption frequency to be greater than the frequency of a prograde 



circular orbit at the ISCO, GMf\^^^/c^ = {7c[x+{c'^Risco/GMY/^]}-^ flBardeen. Press, fc Teukolskv 
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I972I ). then Eqs. f|T6l) - f|T7j) imply that the NS escapes disruption provided that 

7.33km (m/1.4M0) 5 (M/SOMo)^ U + [c'^Risco/ {GM)]^^''^^ x < 0.6894 

1 2 r n 3^ 0-95 

4.59km (m/1.4M0)3 (M/50Mq)3 + [c^Risco/{GM)]2 i x > 0.6894 

(18) 

Reasonable NS models have a maximum radius of ~ 16 km or less, so this criterion will 
be satisfied for a 50Mq IMBH if the spin x < 0.95. For a IOOMq IMBH, the condition is 
satisfied for all spins up to 0.998. As discussed later, we expect IMBHs that grow through 
minor mergers to have only moderate spin x ^ 0.3, so tidal disruption should not occur for 
such IMBHs. 

Although the NS cannot be directly tidally disrupted, tidal oscillations will be excited 
every time the star passes through periapsis. If sufficient energy is deposited into such tides. 



the star could eventually be disrupted through this tidal heating (iFreitad 120031 ). To assess 
whether this effect could be important, we consider the orbital energy lost to leave the star 
on an orbit with periapsis and eccentricity e divided by the binding energy of the star, 
Eorh/Eund- If the inspiral was entirely driven by tidal dissipation, and the tidal energy was 
not efficiently radiated, this would be the ratio of the energy in tidal oscillations to the 
stellar binding energy. Under these assumptions, if this ratio was of the order of 1 or more, 
then tidal heating could disrupt the star. In practice, however, most of the orbital energy is 
lost to gravitational radiation, since, as we see below, tidal oscillations can only be excited 
during the late stages of inspiral. Thus, most of the energy does not go into tidal heating, 
and therefore this ratio would have to be significantly greater than 1 for tidal disruption to 
occur. 



Assuming a Keplerian orbit, this ratio is equal to (lFreitag]|2003l ) 



-1 



where we have assumed that the star has zero kinetic energy at infinity. (Assuming that 
the stellar velocity is 10 km s~^ at infinity changes this result by only 2.3 x 10~^ for a 
I.4M0 NS of radius 10km.) For an inspiral into a Schwarzschild BH, plunge occurs when 
c^rp(l + e) = 2(3 + e)GM; therefore for any eccentricity we have (1 — e)GM/{c^rp) < 1/6 
at plunge. This means that the energy ratio defined in (fT9|) can only be greater than 1 for 
i?* > 12.5 km. Tidal disruption due to heating is very unlikely to occur. This conclusion 
also applies to BHs of moderate spin. For an orbit that is circular at plunge into a BH with 
spin X = 0.35, the ratio -EorbZ-^bind is approximately equal to 1 at ISCO for R^ = 10km. 

If systems existed in which an NS was on a prograde inspiral orbit into a rapidly spinning 
BH, the periapsis at plunge would be much closer to the central body and the energy ratio 
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would exceed unity at plunge. However, the energy ratio would still be small. The radius of 
the ISCO for a BH of spin x = 0-9 is at c^rp = 2.32GM, at which radius i?orb/-£'bind ~ 2 for 
i?* = 10km. The disruption criterion that _Eorb/-Ebind ~ 1 assumes that the orbital energy 
is dissipated entirely by tidal interactions. In practice, the inspiral will mainly be driven by 
GW emission, since most of the orbital energy is lost in the regime where GW emissions 
are quite significant. Tidal dissipation would have to occur on a very short timescale to 
dominate over gravitational radiation reaction effects, and this will not happen in practice. 
We can thus conclude that disruption of the NS due to tidal heating will not occur. This is in 
contrast t o main-sequen ce stars that, being less compact, will be disrupted before reaching 
the ISCO (jFreitag||2003l ). We note that this conclusion does not change when the relativistic 
orbital energy is used in place of the Keplerian expression. 



2.5.2. Tidal Capture 



Although tidal interactions should not shorten the inspiral by causing disruption of the 
NS, if orbital energy and angular momentum of the binary are lost into normal modes of 
the star, the inspiral trajectory will be modified. In principle, this could modify the capture 
rate and the typical eccentricities expected at plunge. Significant oscillations are only likely 
to be excited by tidal interactions if the orbital frequency is comparable to the frequency of 
normal modes in the NS. We can estimate the latter from the frequency associated with the 
free-fall time in the NS: 




5.9UHZ f^Vf^Y 



(20) 



This is just an approximation, but it gi ves the correct order of magnitude for the normal mode 



frequency. iPress fc TeukolskyI (119771 ) computed normal modes using a polytropic stellar 



model with index n = 3. They found an /-mode frequency that agrees with Eq. fl20|) . but 
with a prefactor of 6.2 kHz instead of 5.9 kHz. 

Other stellar modes, in particular gf-modes, can have significantly lower frequency and 
thus will be excited earlier in the inspiral. Press and Teukolsky tabulate frequencies for 
(yf-modes up to (714, which has a frequency a factor of 0.15 smaller than the /-mode. An 
n = 3 polytrope is not a good model for an NS, but it still provides a reasonable estimate 
of the frequency range for thermal (yf-modes. NSs also support crustal ^f-modes that arise 
from chemical stratification a nd co r e q-n iodes that arise from stratification in the number 
densities of charged particles. iFinnI (119871 ) computed frequencies of crustal (/-modes in zero- 
temperature NSs, using a range of stellar models. He fo und that the longest period modes 



had periods of ~ 20ms. iReisenegger fc GoldreichI (119921 ) computed the frequencies of core 
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5'-niodes, and found that these have similar frequencies to the crustal modes. Taking 
as a reasonable maximum for the g-mode period gives a frequency of 20Hz. 



50ms 



Inertial (r-)modes in rotating NSs typical l y hav e frequencies of the order of the spin fre- 
quency of the NS (/ ~ 10 — lOOHz). iHo fc Lail (119991 ) examined the excit ement of r-modes by 
Newto nian tidal driving and found that this was fairly weak. However, [Flanagan fc Racine 
( 120071 ) computed the effect of post-Newtonian gravitomagnetic driving and found that this 
was significantly greater. For rapidly rotating NSs, the inertial-frame frequency can be 
much smaller than th e corotating-frame frequency, which allows /- and p -modes to be ex- 
cited ( IHo fc Lailll999l ). This requires very rapid NS rotation, /rot ~ 500Hz. IHo &: Lail (119991 ) 
examined such modes in the context of comparable mass binaries but concluded that such 
NS spins were unlikely to be found in binary systems. In the IMRI case, where a free NS 
is captured, the NS sp i n cou ld be much higher in prin c iple, m aking these modes potentially 
interesting. iHo fc Lail (119991 ) and iFlanagan fc Racind (120071 ) considered only modes in the 
LIGO frequency range, lOHz < / < lOOOHz, but the mode spectrum extends to lower fre- 
quencies. However, the frequency at which each resonance occurs is a single-valued function 
of the spin of the NS. 



We compare these frequencies to the orbital frequency of a prograde circular orbit at 
radius r: 

M ^- 

+ X 



^^orb = 0.65kHz 



2 

c r 



\GM 



/ M 

VSOMrr 



(21) 



Any NS that comes within a distance ^ 280GM/c^ from an M = 100 M© IMBH will be 
directly captured as a result of GW emission. The additional energy lost in tidal interactions 
could increase this capture cross section. However, for r = 300GM/c^ (cf. § 12.31) . c^orb = 0.13 
Hz, which is much less than the frequency of oscillations of the star. The g-mode frequency is 
2 orders of magnitude higher than the orbital frequency at that radius and so it is unlikely to 
be significantly excited. The g-mode frequencies become comparable to the orbital frequency 
for a Schwarzschild BH when c^r < lOGM. Thus, (^-modes are likely to be excited in the 
late stages of inspiral, but not earlier. As mentioned above, the spectrum of NS r-r aodes and 
the/- and p-mode resonances o f rapidly rotating NSs extend to low frequencies ( IHo fc Lai 
19991 : iFlanagan fc Racind l2007l ). However, the resonant frequencies are determined by the 
NS spin, so it would require extreme fine tuning for a given NS to be captured at precisely 
the periapsis that allows excitement of those modes. The capture rate is unlikely to be 
increased by this mechanism either, although these modes could also be excited later in the 
inspiral. 



Press fc TeukolskyI (119771 ) provide an expression for the energy dissipated in tides in an 



object of mass m and radius R^, that passes a point mass of mass M on a Keplerian orbit 
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with periapsis -Rmin: 

A-E'tidal 




Rn 



2 



(22) 



This expression is integrated over all thermal normal modes, including gf-modes up to qm. 
Once again, this result is based on an n = 3 polytropic stellar model, which is not a good 
model of an NS. However, it should provide an order-of- magnitude estimate for the energy 
lost in thermal modes. In Eq. fl22l) we include only the / = 2 modes, since other modes 
are suppressed by (-R^,/-Rmin)^ <C 1 relative to these modes. We also take the extreme 
mass ratio limit M ^ m. The function T2{r]) behaves as Ty(v) ^ 0.657 7~^'^^ at large t] 



(we have derived this "by eye" from Figure 1 in iPress fc TeukolskyI (119771 )). We can thus 
compute the ratio of the energy dissipated in tides to the energy dissipated in GW emission, 
AEgw = [857rmV(12v^ Mc^)] {GM/R^i^y/^, for an object on a parabolic Keplerian orbit 
with periapsis -Rmin: 



AEtidai f GM R, M Y'-'' f m y''^' 

^■^^[c^Ru^iJ UokmJ UoMeJ [iAmJ 



AE, 



GW 



It is clear that, under these model assumptions, the tidal perturbation to the orbit at capture 
is always much weaker than the perturbation induced by GW emission. For comparison, since 
Aii^GW oc Tp , a 10% increase in the energy lost in a single pass by the central BH increases 
the minimum periapsis required for direct capture by only a factor of 1.1^/^ 1.03 or ~ 3%. 

The above arguments indicate that the excitement of NS modes will neither increase the 
capture rate nor lead to NS disruption during an IMRI. However, orbital energy lost into os- 
cillations could modify the inspiral trajectory by either causing a cumulat ive phase shift in the 



emitted GWs or changing the eccentricity of the orbit in the LIGO band. iFlanagan fc Racine 



(120071 ) calculated the phase difference in the GWs that arises from the excitement of r- 
modes, finding A$ ~ 3.4i?^Q/g/oo^L4"^L4(^i-4 + "^i.4)~"^^^5 where i?io is the NS radius in 
units of 10km, Mi, 4/7711.4 are the masses of the primary /secondary in units of I.4M0, and 
/sioo is the spin frequency (or r-mode frequency) in units of 100 Hz. For an IMRI with 
M = 5OM0, this gives A$ ~ 0.003 if we set Riq = fsioo = mn = 1. Typically we require 
a phase shift of A$ ~ 1 for an effect to be observable, so the excitement of r-modes will 
not leave an imprint on the inspiral. The phase shift induced by the resonant e xcitem ent of 



/- and p-modes in rapidly rotating NSs can be significantly higher. iHo fc Lail (119991 ) quote 
A$ ~ 234m^4^i?fcfm^4/(Mi.4(mi.4 + Mi,^))/-^^^^ for the most extreme case of the (22,2) 
/-mode resonance (with the same notation as before but now denoting the GW frequency in 
units of 100 Hz by fgojioo)- For an M = 50Mq IMRI, this gives A<l> ~ 0.2/^^-^^^. This could 
be a measurable shift if the resonance is excited near 10 Hz. However, the phase shift is only 
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this large for IMBHs at the low-mass end of the IMRI range, and provided that the NS spin 
is tuned to ensure that the resonance is excited near 10 Hz. More work will be needed to 
quantify how large a phase shift would be measurable with Advanced LIGO, accounting for 
correlations between the waveform parameters. 

We can estimate qualitatively what effect tidal dissipation would have on the orbital 
eccentricity and periapsis. The phase-space trajectory that an inspiral follows is determined 
entirely by the ratio dE/dLz- Assuming a Keplerian orbit, we have 



rp (2 v/(l + e)GM - r| dE/dL 
(1 + e) r| dE/dL, - 2(1 - e)y/{l + e)GM 



(24) 



We now suppose that the inspiral was driven entirely by tidal dissipation. Typically the 
dominant excited mode would be an m = 2 mode, for which AL^ = 2/S.E/ujoq, where cjoo 
is the frequency of the mo de (this assumes that the stellar oscillations can be modeled as a 



linear Lagrangian system; iFriedman fc Schutzlll978h . We write 




Woo = 'i/^r' (25) 



where is the radius of the circular (Keplerian) orbit that would have the same frequency 
as the m = 2 mode. With this substitution, equation defines the evolution of r^/ Vc over 
the inspiral. Equations ( 120|) and ( 12T|) indicate that the capture periapsis, r^, will typically be 
much greater than Vc- Solutions with > 2^/^rc are all qualitatively the same, and we show 
a typical example in Figure [21 for capture periapsis of 1000 Vc and a capture eccentricity of 
1. For a IOOMq IMBH, taking ujqq = 6 kHz yields c^r^ ~ 0.5GM, so this figure represents 
a capture at Vp ^ 500 GM/c^, the upper end of the allowed direct capture range for an 
m = 10 Mq BH. The figure shows the inspiral in eccentricity-periapsis space. Under this 
simple model of tidal interactions, the periapsis increases while the eccentricity decreases. 
In practice, the inspiral will be driven by a combination of GW emission and any tidal 
dissipation that occurs. These results suggest that tidal effects would tend to make the 
eccentricities at plunge smaller than they would be for inspirals driven by radiation reaction 
alone. 

Equations (120|) - (!23|) indicate that normal modes are unlikely to be excited during an 
inspiral into an IMBH, although high-order (yf-modes, r-modes, and /-modes in rapidly rotat- 
ing NSs might be excited during the very late stages of inspiral. Thus, we can safely ignore 
the effect of tides on the capture rates. Tidal effects could modify the inspiral, although the 
above calculation indicates that this should not modify our conclusions about the typical 
eccentricities at plunge. The excitement of /-modes might leave a measurable imprint on 
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the GW signal. However, the induced phase shift is only marginally detectable, and this 
mechanism requires the NS to be rapidly rotating. 



3. Event Rates 

In this section we estimate the rate of IMRIs in globular clusters detectable by Advanced 
LIGO. To do this, we must consider three elements: (1) the distance sensitivity of the 
detectors to GWs from IMRIs (and hence the volume of the universe the detectors can see), 
(2) the number density of globular clusters, and (3) the rate of IMRIs per globular cluster. 



3.1. Advanced LIGO IMRI Sensitivity 



For GW sources with known waveforms (or at least waveforms well approximated by 
analytic or nu merical techniques), matched filtering is use d to search for signals in GW 
detector data (jWainstein fc Zubakovlll962l : lAUen et al.ll2005l ). The S/N p of a template h{t) 
in data s{t) collected by a detector that has one-sided noise power spectral density Sn{\f\) 
is given by 

4 r\~sifyh*if)\ 



p 



a 



df, 



(26) 



Snilfl) 

where s{f) is the Fourier transform of the signal s{t), h{f) is the Fourier transform of the 
inspiral template h{t), the asterisk denotes complex conjugation, and a is defined by 



a 



4 



Snilfl) 



df. 



(27) 



This d efinition of S/N follows the normalization of lCutler fc Flanagan! (119941 ) and I Allen et al. 
( I2OO5I ). We place the template h{t) at a canonical source distance of 1 Mpc and choose the 
optimal orientation of the detector to maximize the S/N, and so the maximum distance 
to which a single detector matched filter search is sensitive at a given S/N p is given by 
D = a/p Mpc. (This is the sar ae quantity as the " inspiral horizon distance" used by the 
LIGO and Virgo Collaborations lAbbotLglalJ fl2008h .) 



To compute the sensitivity of a single Advanced LIGO detector to IMRIs, we need 
to compute the quantity cr^ defined in Eq. f l27|) using a particular waveform model. W e 
have done this with waveforms based on BH perturbation theory (IFinn &: Thornd I2OOOI ). 
which are valid in the limit m/M <^ 1. The waveforms, which include non-quadrupolar 
harmonics of the orbital frequency in addition to the dominant quadrupolar harmonic, are 
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described in Appendix \^ where we also discuss the relative S /N contrib uted by the fou r 
lowest harmonics. The noise power spectral density Sn{\f\) was taken from iFritschell (120031 ) . 
GW detectors have an orient at ion- dependent response. The relation between the range R 
(defined as the radius of a sphere whose volume is equal to the volume of the universe in 
which inspiral sources could be detected wi th an S/N threshold of p) and maximum distance 
D at a fixed S/N is given by R = D/2.26 Jpinn fc Chernofjll993l l 



We assume a value of p = 8 for the threshold S/N required for a detection, since this is 
the value typi cally used to comput e Initial LIGO detection ranges for comparable-mass black 
hole binaries (lAbbott et al.ll2008l ). This is a reasonable approximation, as a binary black 
hole inspiral with a total mass of 6 Mq has approximately 500 GW cycles between the 40 Hz 
low-frequency cutoff of Initial LIGO and coalescence — roughly the same number of GW 
cycles that an IMRI signal in Advanced LIGO will have between the Advanced LIGO low- 
frequency cutoff of 10 Hz and coalescence. The threshold will be computed more accurately 
when an IMRI search is implemented and the amount of non-stationarity of the Advanced 
LIGO data is known. If the p = 8 threshold cannot be achieved in practice (or if it can be 
improved), then the detection rates derived below can be scaled appropriately. 

Advanced LIGO will consist of a network of three 4-km detectors. Demanding that 
GWs are found coincident in all three detectors increases the network range by a factor of 
relative to the range of a single detector at a given S /N (due to the lower false alarm rate 
of the network). Fig. [3] shows the range i? of a network of three Advanced LIGO detectors 
for circular-equatorial-orbit IMRIs of m = IAMq objects into a Kerr IMBH of mass M, 
assuming that the network S/N required for a confident detection was p = 8. This is 
equivalent to the range of a single detector with S/N of p = 8/\/3. The x = (non-spinning 
IMBH) line in Fig. |3]is well-approximated by a quadratic fit: 



R ^ ^/tti/Mq 



800 - 540 



M 



100 Mq 



107 



M 



100 Mq 



Mpc. 



(28) 



The scaling of the range in Eq. ( 128|) as -y/m does not follow from the fit, but rather 
from the following reasoning. The amplitude of GWs from IMRIs will scale linearly with the 
mass of the smaller object m, but the number of cycles in the LIGO band will also drop by 
roughly a factor of m. Hence, the total S/N will grow as \/m, so inspirals of more massive 
COs will be seen a factor of ^/m farther away. 



The combination of the spin of the central object and the inclination of the orbital 
plane of the inspiraling particle will have a significant effect on the signal from an IMRI. 
The frequency of the ISCO is much higher for prograde inspirals into rapidly spinning BHs 
than for inspirals into non-spinning holes; the S/N can be strongly enhanced for such orbits. 
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Conver sely, ret r ograde inspirals will have lower S/R. Averaging over random inclination 
angles, iMandell (120071 ) computed the ratio between (1) the detection range for Advanced 
LIGO in a universe uniformly populated by IMBHs of a given mass and spin and (2) the 
detection range in a universe with an equal density of Schwarzschild IMBHs with the same 
mass. He found that the detection range can be enhanced by a factor of 1.7 (3.8) for 
maximally spinning Kerr BHs with M = IOOMq (M = 2OOM0); the increase in the volume 
of observable space and, hence, the event rates, is the cube of these numbers. 

If IMBHs grow mainly by random mergers, they will not be rapidly spinning as the 
contributions of subsequent mergers to the hole's spin largely cancel out. The angular 
momentum imparted to the IMBH by a CO is Lobj oc mM, since the radius at ISCO 
is '"isco oc M. This causes the dimensionless spin parameter of the hole x = Si/M"^ to 
change by ~ Lobj/M^ oc m/M. After ~ M/m such mergers, necessary for the hole to 
grow to mass M, the typical spin of the hole will be ^ ~ \/m/M. More precise calculations 



flHuehes &: Blandfordll2003l : Miller ll2002HMandelll2007h show that the spin of IMBHs involved 
in LIGO IMRIs will rarely exceed x = 0-3 for IMBHs that gained a signifi c ant fr action of 
their mass via minor mergers. For small values of Eq- (24) of iMandell (120071 ) yields a 
correction to the range presented in Eq. (l28l) due to the inclusion of the IMBH spin; the 
detection range in Mpc as a function of M, m, and x is 
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(29) 



This range estimate does not include the cosmological redshift. The redshift due to the 
expansion of the universe decreases the frequency of the GWs. For M ~ 100 Mq IMRIs, the 
redshifted GWs will lie in a less sensitive part of the LIGO noise curve, thereby reducing the 
range. For IMRIs detectable with Advanced LIGO, redshifts are typically < 0.2; for example, 
the inspiral of a 1.4Mq NS into a non-spinning IOOMq IMBH is visible to a redshift of 0.09. 
We estimate that for typical sources, properly including the redshift reduces the Advanced 
LIGO event rate by ~ 10%. 

Advanced LIGO will have several parameters that may be tuned during the operation 
of the detector to optimize the noise power spectral density (PSD) in order to search for 
specific sources. These tunable parameters include the laser power and the detuning phase 
of the signal recycling mirror. If a noise PSD optimized for detections of CO-IMBH binaries 
is used instead of the default PSD assumed in Fig. [3l the range for such sources is increased 
by a factor of ~ 1.5, corresponding to an event rate increase by a factor of ~ 3.5. 
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3.2. Number Density of Globulars with a Suitable IMBH 



The second element in the rate calculation is the number density of globular clusters 
that have an IMBH in the relevant mass range. This is highly uncertain. To contribute 
significantly, a cluster must have had a sufficiently small initial relaxation time to allow the 
formation of an IMBH through some mild runaway process when the cluster was young, 
yet not have formed an IMBH with M > 350 Mq (since this would put IMRls beyond the 
Advanced LIGO frequency r a nge). Recent theo r etical arguments by Trenti and colleagues 
faeggie. Trenti. fc HutI bood : iTrenti et all booj . iTrenti. Heggie. fc HutI bopj . iTrentil booeh 
suggest that dynamically old globulars with large core to half-mass radius ratios have been 
heated by a ~ 1000 Mq IMBH, so these clusters would not contribute to the Advanced LIGO 



IMRI rate. Note, however, that iHurleyl (120071 ) has shown that current observations of the 
core-to-half-light ratios in globulars do not require 1000 BHs in most clusters. Core - 
collapsed globular clusters, which constitute ~20% of all globular clusters (|Phinneylll99ll ). 
may contain IMBHs of the right mass. We will parametrize the unknown fraction of relevant 
globular clusters by some fraction / . Globular clusters have a space density of 8.4 Mpc~^ 
( Portegies Zwart fc McMillanll2000l ). which for h = 0.7 yields 2.9 Mpc~^. Therefore, we will 
use the number density ~ 0.3 (//O.l) Mpc~^. This factor / depends on both the number 
of clusters with an IMRl in the right mass range and the number of clusters that have been 
in a state of core-collapse long enough for the binary hardening mechanism to occur. These 
factors are degenerate, however, since clusters with heavier IMBHs will not be in a state of 
core collapse, as described above. The fraction / also depends on what proportion of the 
objects merging with the IMBH are COs as opposed to main-sequence stars. Our Monte 
Carlo simulations, which were discussed earlier, indicate that this proportion is close to 1. 

The fraction / of globular clusters containing IMBHs may be further lowered by ejections 
of IMBHs from their clusters by recoil kicks imparted to the IMBHs by dynamical processes 
and by gravitational radiation emission. If the kick exceeds ~ 50 km s^^, which is the 
escape velocity from a massive globular cluster, the IMBH will escape from the cluster, 
thereby becoming unavaila ble for future events. Kicks can arise from the process of hardening 



via three-body encounters (IKulkarni. Hut. &: McMillan! Il993l : ISigurdsson fc Hernquistlll993 



Giihekin. Miller, fc Hamiltonlbooi booeh . Giiltekin. Miller and Hamilton (2006) show (cf 



their Fig. 12) that when the seed mass is lOOM©, only about 50% of all BHs grow to 3OOM0 
without being ejected, and this fraction drops to 10% for a seed mass of 5OM0. 

Kicks also arise from GW emission. During the last stages of the merger of unequal mass 



to the resulting BH ( 
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2005 



2007a 



Damour fc Gopakumaiil2006 



Herrmann et al. 2007a: Baker et al. 2006: Gonzalez et al. 



Sopuerta. Yunes. &: Laguna 



20071 ). The most recent results on merger velocity kicks, 
based on numerical relativity, show that the kick velocity for a non-spinning central object 
depends on the symmetric mass ratio 1] = mM/{m + M)^ as \4ick ~ 12000r7^A/l — ^ri{l — 
0.93?7) km s~^ (iGonzalez et al.ll2007al ). The requirement \4ick < 50 km s~^ places an upper 
limit on m of g = m/M < 0.08. 



kick can be a lot higher 
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provide a fit to numerical relativity results that gives the kick as a function of the various 
orbital parameters. This formula indicates that if the IMBH has moderate spin x ^ 0.5 
and the secondary is non-spinning, then we require q < 0.05 to ensure that the IMBH has a 
high probability of remaining in the globular cluster today after undergoing multiple mergers. 
This constraint can be relaxed to g < 0.067 if x ^ 0.3. If the objects merging with the IMBH 
are BHs with a mass of IOMq, this constrains the initial IMBH mass to be M > 15OM0. If 
the merging objects are IAMq NSs, even IMBHs with a seed mass of 5OM0 are safe from 
ejection. 

As argued earlier, mergers with BHs are likely to be important early in the IMBH 
evolution, when its mass is smaller, with NS mergers becoming dominant later. This could 
mean that a significant number of IMBHs were ejected from globular clusters early in their 
evolution. However, without firm knowledge of the initial seed masses of IMBHs or the 
relative number of mergers with BHs and NSs that each IMBH undergoes, it is impossible 
to draw definitive conclusions. We normalize / to 10% in the rate calculations that follow, 
but we emphasize that this quantity is highly uncertain at present. 



3.3. IMRI Rate per Globular Cluster and Event Rate 

The final contribution to the rate estimate is the merger rate per globular cluster. 
Existing numerical simulations of globular clusters suggest that mergers in the subcluster 
of ~ IOMq BHs at the center of the globu lar cluster can lead to the creation of IMBHs 



with masses up to ~ 35OM0 in ~ 10 yr ( lO'Leary et al.ll2006l ). However, the results of 



such simulations are very sensitive to the choice of cluster models and to assumptions about 
kick velocities, the interaction between the BH subcluster and the rest of the cluster, etc. 
Therefore, we present two methods for computing the rate per globular: (1) an upper limit 
independent of cluster model and (2) an estimate based on a more realistic model for cluster 
dynamics. 
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We estimate a theoretical upper limit on the IMR I even t rate in a globular cluster using 



the following method, originally suggested by iPhinneyl (120051 ) . We assume that each globular 



cluster has a BH that grows from M ~ 50Mq to M ~ SSOMq by capturing a sequence of 
COs of identical mass m over the age of the cluster. Then SOOMq/tti captures will happen 
in each globular cluster in ~ 10^° yr. This leads to a rate of (300 MQ)/m x (10^° yr)~^ per 
cluster. 

Although this rate is plausible, it may be a significant over-estimate for several reasons. 
First, it assumes that all the mass that the IMBH acquires in growing from M ~ SOM© 
to M ~ 35OM0 comes from mergers with COs. In practice, the IMBH will also acquire 
mass via gas accretion, and by captures of main-sequence stars and white dwarfs, which will 
be tidally disrupted before becoming significant GW sources but will still add mass to the 
IMBH. Second, this estimate does not include the likelihood that the merger product will 
be kicked out of the cluster through recoil, as discussed in the previous section. Third, this 
estimate assumes that the rate at which the IMBH grows via IMRIs from 50 Mq to 350 Mq 
is constant in time. However, Advanced LIGO can only detect mergers that occurred at 
distances < 1 Gpc, i.e., relatively recently, so the relevant rate is the r ate late in the history 



of the globular cluster, which is likely to be much lower. For example, lO'Leary et al.l (120061 ) 
found in their numerical simulations that the rate dropped from ~ 10~^ to ~ 3 x 10~^° yr~^ 
after 10^° yr for some plausible cluster models. 



For the theoretical upper limit, the total rate is given by ay(M, m, x), where a ~ 
0.3 (//O.l) Mpc"^ (300 M0)/m (10^° yi)-^ is the IMRI rate in the universe, V{M,m,x) = 
(4/3)7ri?'^ is the volume in which Advanced LIGO can see an event, and on overbar, V, 
denotes the average over mass M in the range between 5OM0 and 35OM0. If we take 
/ = 0.1, X = 0-2 as the typical IMBH spin, and all inspiraling objects are 1.4M0 NSs, 
the event rate is ~ 3 yr^^; if / = 0.1, x = 0.2, and inspiraling objects are IOMq BHs, 
the event rate is ~ 10 yr^^. These values are based on the range fit in Eq. (l29l) . so they 
assume that orbital frequency harmonics through m = 4 are included in the data analysis, 
but cosmological redshift and Advanced LIGO optimization are not included. When all of 
these considerations are taken into account, a theoretical upper-limit estimate suggests that 
Advanced LIGO may detect up to 30 IMRIs per year. A similar estimate for Initial LIGO 
shows that because of lower overall sensitivity and a higher low-frequency cutoff (40 Hz for 
Initial LIGO vs. 10 Hz for Advanced LIGO), the upper limit on the Initial LIGO IMRI rate 
is only about 1/1000 events yr~^. 

A more realistic estimate is based on the assumption that the hardening of a CO- 
IMBH binary via three-body interactions represents the primary capture mechanism leading 
to IMRIs. The rate for IMRIs created by this scenario is ~ 3 x 10^^ yr~^ per globular 



- 25 - 



cluster for NS-IMBH IMRIs and ^ 5 x lO'^ yr~i for BH-IMBH IMRIs [see § |2ll]. Hence, 
the NS-IMBH IMRI rate in the local universe is a ^ 10"^ (//O.l) Mpc"^ yr~\ while the 
BH-IMBH IMRI rate is a ^ 1.5 x lO'^ (//O.l) Mpc"^ yr^^. If we assume that all IMBHs 
have a mass ~ lOOM© and / = 0.1, this yields an Advanced LIGO rate of one IMRI per 3 
years if the typical CO is an NS or 10 IMRIs per year if the typical CO is an m = 10 Mq BH. 
If the interferometer is optimized for the detection of IMRIs, the NS-IMBH and BH-IMBH 
rates are increased to 1 and 30 events yr~^, respectively. 

In addition to detections of inspirals. Advanced LIGO could also detect the ringdown 
of an IMBH following a merger. This possibility is discussed in Appendix [Bl 



4. Effect of Eccentricity on Matched Filter Searches 

As discussed in § 13.11 matched filtering is used to search for GWs with known waveforms 
in detector noise. In order to be an optimal search technique, the matched filter requires 



accurate templates that correctly model the signals being sought (IWainstein &: Zubakov 



19621 ). Since source parameters (e.g., the masses and the IMBH spin) can vary, the matched 
filter is constructed for a "bank" of templates: a set of waveform models that depend on the 
parameters that character ize the source. Th e accuracy of a template bank is characterized 



by the fitting factor (FF) (lApostolatoslll996l ). which measures the overlap between the GW 
signal and the nearest template. A fitting factor close to unity indicates that the templates 
are accurate for detection of the desired signals. A fitting factor less than unity will mean 
that we are unable to detect a fraction (1 — FF^) of the theoretically detectable events. (The 
quantity 1 — FF is often referred to as the mismatch.) To search for signals, template banks 
are constructed so that the mismatch between any desired signal and the nearest template 
does not cause an unacceptable loss in S/N (typically FF ^ 0.97 for LIGO). 

In this section we examine the effect of eccentricity on searching for IMRI signals in 
Advanced LI GO detectors. Th e effec t of eccentricity on the fitting factor was previously 



examined by iMartel fc Poissod (119991 ). and it was found that the fitting factor between a 
circular and eccentric waveform template was high provided that e < 0.2. However, their 
results do not apply directly to IMRIs since they computed fitting factors only for binaries 
with mass ratios close to 1, and used the first-generation LIGO noise curve. 

We consider a matched- filter search for IMRIs and determine the loss in S/N (and hence 
range) if eccentricity is not included in the template bank; i.e., circular templates are used 
to search for potentially eccentric waveforms. We compute the fitting factor as follows. The 
template h{t) appearing in the expression for the matched filter S/N p [Eq. fl2B|) ] depends 
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on a number of parameters characterizing the source, such as the masses of the binary and 
the time of arrival of the signal. We denote these parameters A and define the ambiguity 
function ^(A) by 

^(X) = , <--l'"^"» . (30) 



(s\s){h{\m\)} 

where {a\b) is the matched filter inner product given by 



We separate the parameters A into A = (to,0o,^), where to and 0o are the time of 
arrival and phase of the binary, respectively. In the case of circular equatorial binaries, it is 
trivial to maximize over the parameters to and 0o analytically (the phase by projecting the 
signal onto two orthogonal basis vectors and the time by a Fourier transform), and so these 
are called "extrinsic parameters." The remaining template parameters 6, which include the 
binary masses, eccentricity, and IMBH spin, determine the shape of the waveform and are 
known as "intrinsic parameters." For circular inspiral templates, the ambiguity function A 
reduces to the overlap O, given by 

0(^-') = max^Jfil^, (32) 

^{s\s){h{e)\h{e)) 

The fitting factor is given by the maximum of the overlap function over the remaining 
parameters 

FF = maxC»(^). (33) 

e 

For the signal s{t) and template h{t) we use numerical kludge waveforms. This is a family 
of waveforms that were constructed as models for extreme mass ratio inspiral systems, in 
which m/M ^ 1. The waveform family is constructed by first computing an accurate phase- 
space trajectory by integrating prescriptions for the evolution of the orbital elements (the 
orbital energy, angular mome ntum, and Carter constant, or equivalently the orbital radius. 



eccentricity, and inclination) (IGair fc Glampedakid 120061 ). The orbit of the small body is 



then calculated by integration of the Kerr geodesic equations along the sequence of geodesies 
defined by the phase-space trajectory. Finally, a kludge waveform is generated from the 
orbit by applying weak- field emission formulae ( Babak et al. 2007h. This waveform family 



predicts the inspiral rates for nearly-circular orbits very well (jCair &: Glampedakis 2006) 



and has been shown to be extremely faithful (overlaps in exc ess of ~ 95% over much of the 



parameter space) to more accurate perturbative waveforms (IBabak et al.l 120071 ). Although 
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the mass ratio of an IMRI system is probably too high to make these waveforms accurate as 
search templates, they should provide reliable predictions of the fitting factor. 

For these calculations we used M = 100 Mq for the IMBH mass, m = 1.4 Mq for 
the companion mass, and considered two spin values x = and y = 0. 2. We used the 



Advanced LIGO power spectral density 5*^(1/1) given by iFritschell (120031 ). As discussed 
above, to compute the fitting factor one must maximize over the parameters 6 of the template. 
However, we find that even without maximizing over the intrinsic parameters, the overlap 
(and hence the fitting factor) between circular and eccentric templates is greater than 0.99 
for eccentricities e < 0.01, i.e., for more than two-thirds of IMRIs formed by direct capture 
(the mechanism likely to give the largest eccentricities). Since we expect that most of the 
IMRI systems will have eccentricities significantly less than e = 0.01 by the time they have 
entered the Advanced LIGO band, eccentricity will be negligible for data analysis and circular 
templates may be used to search for these systems. 

Fig. m shows the overlap between eccentric signals and circular templates for prograde 
equatorial inspirals and eccentricities greater than 0.01. Analysis of inclined inspirals demon- 
strates that the overlaps between eccentric signals and circular templates remain greater than 
0.99 for eccentricities e < 0.01 and greater than 0.93 for eccentricities e < 0.05. Although the 
overlap decreases for eccentricities greater than 0.01, we anticipate higher values of the fitting 
factor when we maximize over the other intrinsic parameters. An interesting question will 
be to determine whether eccentricities greater than 0.01 can be measured (and thus be used 
to investigate the relative prevalence of the various capture mechanisms) or if eccentricity is 
degenerate with masses and the other intrinsic parameters. 



5. Summary 

In this paper we have discussed a potential source of GWs for ground-based interferom- 
eters: the intermediate mass ratio inspiral of a stellar mass CO (an NS or BH) into an IMBH 
in the center of a globular cluster. For IMBHs with masses in the range 50 — 35OM0, the 
GWs emitted will be at frequencies in the Advanced LIGO band. We have shown that Ad- 
vanced LIGO should be able to detect the inspiral of a IAMq NS into an IMBH at distances 
up to 700 Mpc, depending on the mass and spin of the IMBH. Assuming that all IMBHs 
were grown by CO-IMBH mergers gives an upper limit on the Advanced LIGO event rate 
of ~ 10 yr^^. We have shown that if the inspiraling CO is an NS, a more likely estimate 
of the rate is one event per 3 years, while the rate for BH-IMBH IMRIs could reach the 
upper limit. If Advanced LIGO is optimized for detections at low frequencies, the event rate 
estimates would increase by a factor of ~ 3.5. 



- 28 - 



We have also discussed four mechanisms by which such IMRI systems could form: (1) 
binary hardening via three-body interactions, (2) hardening via Kozai resonance, (3) direct 
capture, and (4) tidal capture of a main-sequence star. In all four cases, we find that the 
residual eccentricity when the inspiral enters the LIGO sensitivity band will be small. Finally, 
we have estimated the sensitivity of Advanced LIGO to the eccentricity of IMRI systems. 
We have found that the eccentricities we expect are negligible for data analysis, and therefore 
circular-orbit templates may be used to search for IMRI binaries in Advanced LIGO. 

IMRIs are a somewhat speculative source of GWs, since evidence for the existence of 
IMBHs is not yet conclusive. The body of evidence is steadily growing, however. Since little 
is known about the abundance of IMBHs in the universe, the event rates presented here 
are naturally somewhat uncertain. However, our results are sufficiently promising to make 
IMRIs a source worth searching for in Advanced LIGO data. If IMRI events are detected 
with Advanced LIGO, these will provide irrefutable evidence for the existence of BHs with 
intermediate mass and will provide information on the mass and spin of IMBHs, plus the 
eccentricities of the inspiraling objects. This information will be very useful for constraining 
models of IMBH formation and growth and for exploring stellar dynamics in the centers of 
globular clusters. 
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APPENDIX 



A. Waveforms and Signal-to-Noise Ratio Calculation 

To compute the range to which a source can be seen, as presented in § 13.11 we must 
evaluate the S/Ns of typical sources. To do this requires a model of the waveform. In the 
weak field, waveforms may be well approximated by post-Newtonian results. The leading- 
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order post-Newtonian result takes the system to be a Keplerian binary and estimates the 



gravitationa l radiation from the leading-order quadrupole formula (jPeters fc Mathewslll963 



Peters! 1 1 9641 ) ■ This predicts h(f) oc /~^''^0(/ — /isco); where the step function 9 is included 
to ensure that the radiation cuts off at /isco; the GW frequency at the innermost stable 
circular orbit of the binary. The post-Newtonian results are a weak field expansion and are 
only valid where velocities are much less than the speed of light. As a consequence, the 
leading-order post-Newtonian waveforms over-predict the S/N of an IMRI source, since they 
effectively spend too many cycles at each frequency as the ISCO is approached. 

An alternative GW model can be obtained from perturbation theory, by expanding 
in terms of the mass ratio, m/M, assumed to be small. The IMRI systems considered in 
this paper lie somewhere between these two extremes — the mass ratio is not quite small 
enough to use perturbative techniques, but the source spends a long time in the regime 
where post-Newtonian results are not valid. Waveform models have not yet been developed 
specifically for IMRI systems. However, by the time Advanced LIGO comes online, it is 
likely that models will have been constructed by co mbining post-Newton i an an d perturbative 
techniques. This is discussed in more detail by lAmaro-Seoane et al.l (120071 ). If accurate 
waveforms are not available, we will require sources to have higher S/Ns to be detected, thus 
reducing the ranges from the values that we quote. However, th e loss in S/N from using a n 
inaccurate template is likely to be only a few tens of percent (lAmaro-Seoane et al.l 120071 ). 
which is considerably smaller than the uncertainties in the astrophysical mechanisms that 
govern the event rates we are computing. 

Out of the set of currently available waveform families, we believe that the most accurate 
S/Ns will come from the perturbative waveforms. Although the perturbative waveform will 
not be a precise model of the true waveform, the total energy content of the GWs will be 
roughly correct since the perturbative methods use a reliable model of the spacetime close to 
the central BH. To generate the range estimates quoted in this paper, we therefore computed 
the S/N via a perturbative model, as described below. 



Finn Sz Thornd (120001 ) used perturbation theory to compute the S/N, averaged over 



sky location and source orientation, contributed by the lowest four harmonics of the orbital 
frequency for circular, equatorial inspirals into Kerr BHs. Their calculation is accurate in the 
sense that it is based on perturbation theory, but it relies on three assumptions: (1) the orbit 
is in the extreme mass ratio limit, i.e., m/M <^ 1; (2) the orbit of the small body is circular; 
and (3) the orbit of the small body is equatorial. Assumption (2) is valid for our case, and 
assumption (1) is probably sufficiently accurate (the mass ratio here is intermediate while 
not extreme). Assumption (3) is not necessarily va lid, but we can derive results for both 
prograde and retrograde equatorial orbits from the iFinn &: Thornd (120001 ) waveforms and 
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then average over possible orbital inclinations of the inspiraling object by assuming that the 
effect of averaging is the same as it is for the leading-order post-Newtonian model (IMandel 
2007h . 



It is conventional to use m to denote harmonic number when discussing harmonics of 
the azimuthal frequency. However, in this paper we use k to avoid confusion with the mass of 
the CO. The S/N contributed by the kth harmonic of the orbital fr equency, A = kuJnrh/C^ ii), 



and averaged over sky location and source orientation is given by (IFinn fc Thornell2000l ) 



pI 



[hcAfk)y 

fkSnifk) 



dlnfk 



(Al) 



where Sn{f) is the one-sided power spectral density of the detector noise and hc^fk) is the 
characteristic amplitude of the fcth harmonic when it passes through frequency fk- This 
reduces to the earlier expression fl27|) via the substitution 2h{f) = J2k hc,k{f)/ f ■ The char- 
acteristic amplitude is related to the energy radiated to infinity in each harmonic and is 
given by 
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for m > 2. (A3) 



Here Vl = GMuorb/c^ is the dimensionless orbital angular frequency and tq is the distance 
to the source. The relativistic correction. Tick, can be written as 



n 



c,k 



ook- 



(A4) 



In this expression. A/" is the relativistic correction to the number of cycles spent near a 
particular frequency, and Sook is the relativistic correction to the rate of energy lost to 
infinity in harmonic m. These corrections can be computed via integratio n of the Teukolsky- 
Sasaki-Nakamura equations and are tabulated in iFinn fc Thornd (120001 ). We note that the 
various corrections are defined relative to their Newtonian values. 



Using the results of lFinn fc Thornd (l2000l ). we can compute the total sky- and orientation- 
averaged S/N ptot contributed by the lowest four harmonics of the orbital frequency from 
the time the source enters the detector band (when = lOHz) until plunge, for various 
spins and masses of the central BH: 



Ptot 



HA- 

k=l 



(A5) 
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This S/N was used to derive the range formulae presented in § 13. 1[ We can also compute 
the leading-order post-Newtonian S/N by including only the quadrupolar k = 2 mode and 
setting the correction Hc,2 = 1- We find that for x ^ 0.5 and 50Mq < M < 25OM0, the 
p ost-Newtonian S/N is t ypically an overestimate by a factor of ~ 1.4. We note that the data 



m 



Finn fc Thornd (120001 ) do not extend to the full range of radii needed for these calculations. 
Where necessary, we extrapolated their results to larger radius using appropriate power laws. 
We have verified that the results are insensitive to the exact form of this extrapolation. 

The simplest template to use to detect a circular inspiral would include only the dom- 
inant, quadrupolar, component of the orbital frequency. It is useful to estimate how much 
S/N we would lose by ignoring higher harmonics. For circular inspirals in the equatorial 
plane of a Kerr BH, the fraction of the total energy radiated during an inspiral from infinity 
that is radiated between a certain Boyer-Lindquist radius and plunge effectively depends 
only on the ratio of the initial radius to the radius of the innermost stable circular orbit, 
rj/risco(x), and is otherwise independent of x- Here x is the central BH spin as usual, and 
Tisco is the radius of the innermost stable circular orbit, as given in Eq. (fT7|) . The energy 
radiated in higher harmonics of the orbital frequency is suppressed relative to that in the 
dominant k = 2 harmonic by powers of M/r. As the BH spin increases, risco/M 1 for pro- 
grade orbits, and so a larger fraction of the energy is radiated in the regime where r ~ M. 
We would therefore expect higher harmonics to contribute most significantly to the total 
energy fiux for prograde inspirals into BHs with large spins. We computed the fraction of 
the total energy radiated into each harmonic as a function of the BH spin, w hile the particle 



i nspir als from r = 10 Tisco to r = Tisco- This is the range of radii for which iFinn fc Thorne 



(120001 ) tabulate data, and in this range ~ 85% of the total energy is radiated in any circular 
equatorial inspiral. The energy fractions are shown in Figure [51 We see that for |x| < 0.3, 
which is the expected IMBH spin range if the IMBH grows via minor mergers, ~ 8% of the 
energy is radiated into harmonics other than the dominant k = 2 harmonic, and most of this 
energy goes into the k = 3 harmonic. 

The contribution of a harmonic to the S/N of a source depends not only on the energy 
that goes into that harmonic, but also on the shape of the noise curve — higher harmonics 
enter the detector band earlier, contribute their signal at frequencies where the noise power 
spectral density is lower, and therefore have an enhanced contribution to the S/N. Figure [HI 
shows the relative S/N contributed by each harmonic, defined as pk/p tot) as a function of 
IMBH mass, for various IMBH spins. Note that this result does not depend on the mass 
of the inspiraling CO, since we are working in the extreme mass ratio limit. We see that 
for prograde inspirals, we can lose ~ 10 — 25% of the S/N by using templates containing 
the k = 2 mode only, but this is mostly recovered by including the k = 3 mode in the 
search templates. (We can lose up to ~ 50% of the S/N by using only k = 2 templates 
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for retrograde inspirals into high- mass IMBHs, but the S/Ns for such events are very small, 
making their detection unlikely.) 

The S/Ns computed from these perturbative waveforms are not totally accurate for the 
reasons given earlier. Corrections will include finite-mass effects, contributions from the spin 
of the small BH and the effect of A; > 4 harmonics of the orbital frequency. It is clear from 
Figure [5] that for larger spins, a significant amount of energy goes into harmonics with k > 4. 
These harmonics spend even longer in band and so their inclusio n would increase the S /N. 
However, we cannot compute their contribution to the S/N since iFinn fc Thorne J2nnnl ) do 



not tabulate these contributions separately. Overall, the S/Ns computed here should be 
accurate to ~ 10% and will be more accurate than those computed from the leading-order 
post-Newtonian waveforms. 



B. Ringdowns 

Following the coalescence of an IMBH with a CO, the BH enters the ringdown phase, 
characterized by oscillations of its quasi-normal modes, particularly the dominant I = 
m = 2 m ode. For IMRIs, the tota l energy emitted in GWs during the ringdown is ~ 



O.SmVM (iFlanagan fc Hugheslll998l ). which is a factor of 0{m/M) sma ller than the tot al 



energy emitted over the inspiral. However, the ringdown GW frequency (lEcheverrial Il988l ) . 



is higher than the ISCO frequency and is therefore closer to the minimum of the Advanced 
LIGO noise power spectral density for the typical masses under consideration. For this 
reason, ringdowns may be detectable by Advanced LIGO despite their lower energy content. 
This is particularly true if m > 10 Mq BHs, rather than NSs, are common as inspiraling 
companions, since the range for ringdowns scales as at low redshifts. Moreover, ringdowns 
will be the only way to detect CO coalescences with slowly-spinning IMBHs with masses 
above 350 Mq, since inspirals into such massive IMBHs will produce GWs at frequencies 
below the detector low- frequency limit. 

The typical Advanced LIGO ringdown- wave ranges (in terms of luminosity distance) as 
a function of IMBH mass are plotted in Fig. [7] for several choices of inspiraling object mass 
and IMBH spin. Because some ranges reach out to significant redshifts (up to 2; ~ 0.5), 
the effect of redshifting is already included in these ranges, unlike in Fig. [31 Redshifting 
also explains why the range does not scale strictly as m^, as high redshifts bring the GW 
frequency at the detector down into the region where the interferometer is less sensitive. 
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The astrophysical rate of ringdowns per cluster is greater than or equal to the rate of 
IMRIs, since every IMRI culminates in a merger and ringdown (but ringdowns could follow 
coalescences without observable inspirals, i.e., those with direct plunges). The distance 
sensitivity to ringdowns following inspirals of 1.4 Mq NSs is probably too low to make them 
detectable by Advanced LIGO: the total detectable event rate for NS-IMBH ringdowns is 
~ 20 times lower than the event rate for NS-IMBH inspirals if the IMBH mass is M = 100 Mq 
and spin is x = 0.3. However, Advanced LIGO will be considerably more sensitive to 
ringdowns than to inspirals in other mass ranges. For example, ringdowns from 10 Mq + 
300 Mq coalescences could be detected in a volume ~ 200 times greater than the detection 
volume for inspirals from these coalescences; if all IMBHs had mass M = 3OOM0, and all 
COs were m = 10 Mq BHs with coalescence rate equal to ^ 5 x 10~^ yr~^ per cluster as in 
§ 12.11 then the total detectable IMRI ringdown event rate would reach ~ 20 yr~^ . Thus, if 
our expectations about the likely masses involved in IMRIs are incorrect, and coalescences 
of COs with higher masses with more massive IMBHs are common, searches for ringdown 
waves can provide a useful back-up to IMRI searches. 
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Fig. 1. — Eccentricity at /gw = 10 Hz as a function of the periapsis at capture, for a CO 
inspiraling into an IMBH of mass M = 100 Mq. The eccentricity at capture is set to 1, and 
the eccentricity at Vp ^ 16 GM/c^, where /gw = 10 Hz, follows from Eq. (Q. 
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Fig. 2. — Tidal-dissipation-driven inspiral in phase space for an inspiraling star with initial 
eccentricity of e = 1 and initial periapsis = 1000 r^. The plot shows eccentricity on the 
horizontal axis and the ratio r^/r^ on the vertical axis. The radius characterizes the 
frequency of normal modes in the star as defined by Eq. (l25ll . 
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Fig. 3. — Range of a network of three Advanced LIGO detectors for the circular-equatorial- 
orbit inspiral of a 1.4 Mq object into an IMBH, as a function of IMBH mass M. The three 
lines show IMRI spins of x = 0-2 (dashed), (solid), and —0.2 (dot-dashed). Positive x 
means prograde orbit; negative x means retrograde. The quadratic fit given in Eq. (128|) is a 
fit to the X = line. 
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Fig. 4. — Overlaps O between a circular template h(t) and signals s(t) with varying eccen- 
tricities, e. For both signal and template, the intrinsic parameters 6 = {M = 100 MQ,m = 
1.4 Mq, X, e) are kept constant, with maximization performed only over time of arrival and 
phase. The overlaps for two values of x a-re shown. 
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Fig. 5. — Fraction of the total energy radiated into each harmonic of the orbital frequency as 
the particle inspirals in a circular equatorial orbit from lOrisco to risco- This energy fraction 
is shown as a function of BH spin. 
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Fig. 6. — Ratio of sky- and orientation-averaged S/N contributed by the lowest four har- 
monics of the orbital frequency to the total sky- and orientation-averaged S/N contributed 
by the lowest four harmonics, as a function of the central BH mass, for circular equatorial 
orbits. The harmonics arc indicated by different line styles — k = 1 (dashed), k — 2 (solid), 
k = 3 (dotted), and k = A (dot-dashed). Lines are shown for three different BH spins, X = 0, 
X = 0.5 and x = —0.5 (i.e., retrograde inspirals into a x = 0.5 BH), indicated by crosses, 
circles, and plus signs, respectively 
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Fig. 7. — Range of a network of three Advanced LIGO detectors for the ringdown of an 
IMBH following a merger with a CO. The luminosity-distance range in Mpc is plotted 
as a function of IMBH mass M; cosmological redshift is included. Dashed lines denote 
m — 1.4 Mq inspiraling NSs, with plus signs corresponding to IMBH spin x — 0-3 and 
crosses to x = 1. Solid lines denote m = 10 Mq inspiraling BHs, with circles, squares, and 
triangles corresponding to spins X = 0) X = 0-3, and % = 1, respectively. The dotted line 
with asterisks denotes m = 20 Mq BHs spiraling into an IMBH with spin % = 0.3. 



